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ABSTRACT 

Context. The O-type object 9 Sgr is a well-known synchrotron radio emitter. This feature is usually attributed to colliding-wind binary 
systems, but 9 Sgr was long considered a single star. 

Aims. We have conducted a long-term spectroscopic monitoring of this star to investigate its multiplicity and search for evidence for 
wind-wind interactions. 

Methods. Radial velocities are determined and analysed using various period search methods. Spectral disentangling is applied to 
separate the spectra of the components of the binary system. 

Results. We derive the first ever orbital solution of 9 Sgr. The system is found to consist of an 03.5 V((f + )) primary and an 05- 
5.5 V((f)) secondary moving around each other on a highly eccentric (e = 0.7), 8.6 year orbit. The spectra reveal no variable emission 
lines that could be formed in the wind interaction zone in agreement with the expected properties of the interaction in such a wide 
system. 

Conclusions. Our results provide further support to the paradigm of synchrotron radio emission from early-type stars being a mani- 
festation of interacting winds in a binary system. 
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1. Introduction 

It has been known for a while that a subset of the early-type stars 
of spectral-type O or Wolf-Rayet display a non-thermal radio 
emission in addition to the thermal free-free emission produced 
in their stellar winds (Abbott et al. I1984I for a recent review 
see Blomme I201 11 1. Dedicated multi-wavelength investigations 
of these objects have shown that they are mostly binary systems, 
hence suggesting that the acceleration of the relativistic electrons 
responsible for the synchrotron radio emission is related to the 
existence of hydrodynamical shocks between the stellar winds 
of the components of a massive binary system (e.g. Dougherty 
& Willia ms [20001 De Becker |2007] Naze et al. 1201 Oal Sana et 
al. l20TlT) . 

One of the first O-type stars discovered to display such a 
synchrotron radio emission was 9 Sgr (= HD 164794, Abbott et 
al. [1984). This star is located inside the very young open clus- 
ter NGC 6530 and contributes significantly to the ionization of 
the Lagoon Nebula (M 8). Until recently, 9 Sgr was considered 
as presumably single. However, its X-ray properties as well as 
episodic variations of some optical line profiles over long time- 
scales (Rauw et al. !2002al i suggested that it could be a binary 
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system. In 2004, we observed another episode of optical line 
profile variations that we then interpreted as the partial deblend- 
ing of the spectra around quadrature phase in a highly eccentric 
binary system with a period of about 8 - 9 years (Rauw et al. 
2005Q 

Synchrotron radio emission requires the presence of rela- 
tivistic electrons and a magnetic field. Concerning the latter, 
Hubrig et al. (2008 ) claimed the detection of a magnetic field 
of about 200 G in 9 Sgr, although this field was only detected 
above the 3-cr confidence level in one out of three observations 
of this star and the other two observations yielded a longitudinal 
field of opposite sign (which could be due to a different rota- 
tional phase). As to the origin of the relativistic electrons, the 
present paper aims at establishing the binary properties of 9 Sgr 
to pave the way for future modelling efforts (e.g. Blomme et al. 
120101) . 



2. Observations 

Since our discovery of the probable binarity of 9 Sgr, we moni- 
tored this star with a variety of telescopes and instruments. The 



1 There exists some confusion in the literature about the orbital pe- 
riod of 9 Sgr. Indeed, the statement by Hubrig et al. J201 lal 1201 lbb 
of a 2.4 yr orbital period, misquoting respectively Naze et al. (2010a) 
and Naze et al. (2010b), is definitely wrong. Actually, whilst the work 
of Naze et al. (20103 refers to another non-thermal radio emitter, 
CygOB2 #9, which has indeed an orbital period of 2.4 yr, this result 
as well as the paper of Naze et al. ( 2010b I are totally unrelated to 9 Sgr. 
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Fig. 1. Various lines in the spectrum of 9 Sgr observed at two different epochs. The lower spectrum was taken on 10 July 2004 
(HJD 2453196.586) whilst the uppermost spectrum was obtained on 5 October 2006 (HJD 2454013.526). From left to right, these 
lines illustrate the various behaviours that exist in the spectrum: lines predominantly formed in the atmosphere of the primary (N v 
AA 4603, 4620), lines not showing any clear binary signature (He n A 5412), lines belonging to the secondary (C m A 5696) and lines 
displaying an SB2 signature (C iv AA 5801, 5812). In each panel, the dashed red lines indicate the rest wavelengths of the lines. 



dataset analyzed in this paper (see Table A. 1 1 was taken over a 
timespan of twelve years. 

The vast majority of our data was obtained at ESO at the 
La Silla observatory (Chile) with the FEROS echelle spectro- 
graph on the ESO 1.5 m telescope until October 2002 and on 
the MPG /ESO 2.2 m telescope afterwards. FEROS (Kaufer et 
al. [1999) has a spectral resolving power of 48000 and covers 
the spectral domain from slightly below 3800 to 9200 A. The 
data were processed with an improved version of the FEROS 
pipeline working under MIDAS. The merged spectra were nor- 
malized over the spectral domain from 3900 to 6700 A by fitting 
a spline function through a large number of carefully selected 
continuum windows. 

Fifteen additional spectra were obtained with the Coralie 
spectrograph (R = 55000) mounted on the Swiss 1.2 m Leonhard 
Euler Telescope at La Silla. Coralie is an improved version of the 
Elodie spectrograph (Baranne et al. 119961) and covers the spec- 
tral range from 3850 to 6890 A. The data were first reduced with 
the Coralie pipeline and subsequently normalized in the same 
manner as the FEROS data. 

One spectrum was obtained with the ESO-New Technology 
Telescope at La Silla equipped with the EMMI instrument in 
echelle mode (Dekker et al. 1986). EMMI was operated with 
grating #9 and grism #3 as a cross-disperser to provide a resolv- 
ing power of R = 7700 between about 4000 and 7200 A. The 
data were reduced with the echelle context of the MIDAS soft- 
ware and the normalization was done order by order using 3rd 
degree polynomials. 

Several observations of 9 Sgr were taken with the UVES 
echelle spectrograph (D'Odorico et al. 2000 ) on UT2 at the ESO 
Cerro Paranal observatory (Chile). UVES was used with the 
standard settings at 4370 A (blue arm) and either of the standard 
settings at 5800, 6000 or 7600 A (red arm). The slit width was 
set to 1 arcsec, yielding a resolving power around 40000. The 



data were reduced with the UVES CPL pipeline and normalized 
in the same way as the FEROS spectra. 

Another set of spectra was obtained with the Espresso 
spectrograph at the 2.12 m telescope at the Observatorio 
Astronomico Nacional of San Pedro Martir (Mexico). The re- 
solving power of this instrument is about 18000. The data re- 
duction was done with the echelle context of the MIDAS soft- 
ware. Due to the strongly peaked blaze, the normalization was 
performed on a limited number of orders (around the most im- 
portant spectral features) using spline functions adjusted to care- 
fully chosen continuum windows. 

We finally used a set of archive spectra collected with the 
REOSC-SEL echelle spectrograprj^Jon the 2.15 m Jorge Sahade 
telescope of Complejo Astronomico El Leoncito (Argentina). 
Again, normalization was found to be challenging due to the 
strongly peaked blaze function and was thus restricted to spe- 
cific orders. 



3. Data analysis 

Figure[T] illustrates some selected parts of the spectrum of 9 Sgr 
at two different epochs. The changing position of the lines, and 
in some cases, the deblending of the lines are clearly seen. 

As a first step, we thus measured the radial velocities (here- 
after RVs) of individual lines by fitting Gaussian profiles. For 
those lines that belong to only one component of the binary 
system (see below), a single Gaussian profile was fitted. For 
those lines that display a clear SB2 signature, we fitted two 
Gaussians whenever the lines were sufficiently deblended to pro- 
vide a meaningful fit. This allowed us to obtain the RVs to be 
used in the orbital solution (see Se cts.|3.1| and |3.2| > as well as in 
the disentangling process (see Sect.|3.3|> 
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We note a variety of behaviours (see Fig.[TJ. Apart from 
some changes in the line width and depth, no clear signature 
of the orbital motion is found for the He n AA 4200, 4542, 4686, 
5412 lines. The Hei AA4471, 5876, Om ,15592, Civ /U5801, 
5812 absorption lines show a clear SB2 signature in our data 
taken in 2004 near maximum radial velocity separation, with the 
blue-shifted component being the stronger one (we shall refer 
to this star as the secondary in the following. The same holds 
for a weak absorption at 4212 A that we tentatively associate 
with Si iv. On the other hand, the He i A 47 1 3 absorption and C m 
A 5696 emission move in phase with the secondary star but show 
no obvious SB2 signature. Finally, the Si iv A 41 16 emission, the 
Nv^ 4603, 4620 absorptions as well as the N iv AA 5200, 5205 
absorptions display an SB 1 signature in phase with the primary. 
Therefore, we compute the RVs of the secondary as the mean of 
the measurements of the secondary components in He i AA 4471, 
5876, Om ,15592, Cm ,15696 and Crv ,1,15801, 5812. All of 
these lines, except Cm A 5696, display an SB2 signature near 
maximum separation, but the secondary's line always dominates 
the blend in such a way that we are rather confident that the mean 
value of these measurements should be a reliable indicator of the 
secondary's RV. In the same way, the primary RVs are obtained 
from the mean of the RVs of the Si iv A 41 16 and Nv,M4603, 
4620 lines, which are probably free of blending with the sec- 
ondary's lines. The results are listed in Table A. 1 To quantify 



the uncertainty on the RV measurements, we have evaluated the 
1-cr dispersion about these means for each star, after correcting 
for systematic shifts between the radial velocities of some lines 
(see also Sect. |3.2] >. 

The same partial line splitting as seen in our 2004 data (see 
Fig.[T]i was observed in 1987 by Fullerton (1 19901 1 and is also 
present in the REOSC spectra taken in 1995. This suggests an 
orbital period around 8.5 years, making 9 Sgr one of the spectro- 
scopic O-type binaries with the longest orbital period so far. 

Williams et al. ( 1201 11 1 recently presented a set of 16 RV mea- 
surements of 9 Sgr, obtained over 13 consecutive nights in May- 
June 2004, i.e. near maximum RV separation. These observa- 
tions were taken in the blue (4058 - 4732 A) with the Ritchey- 
Chretien spectrograph at the 1.5 m telescope at CTIO with a re- 
solving power of 2750. This set-up lacks the resolution that is 
needed to observe the SB2 signature even at maximum RV sep- 
aration, explaining why Williams et al. J201 11 1 did not observe 
the SB2 signature. 

3.1. The orbital period 

To constrain the orbital period of the system, we used the gen- 
eralized Fourier periodogram technique proposed by Heck et al. 
( 1985 ) and Gosset et al. (2001). This method was applied to var- 
ious combinations of the secondary RVs (our data, our data com- 
bined with the Fullerton 1990 data points) as well as the differ- 
ence between the secondary and primary RVs. The best period is 
found to be 3165 ± 1 10 days. The RV data folded with this period 
reveal a highly non-sinusoidal radial velocity curve indicating a 
rather large eccentricity. To check this result, we also applied the 
trial period method of Lafler & Kinman ( |19651 l. This method is 
a priori less sensitive to the shape of the RV curve. Again, using 
various combinations of the RVs, we obtain a best estimate of 
the period of 3181 ± 1 lOdays, which agrees quite well with the 



value obtained from the Fourier technique, given the uncertainty 
of the periocQ These estimates of the orbital period were subse- 
quently used as input to the orbital solution code. We stress that 
the REOSC and Fullerton (1990 ) data actually have very limited 
weight in the period determination process. Indeed, not includ- 
ing these data changes the best estimate of the period by only 
about 15 days. However, including the old data points actually 
reduces the uncertainty on the period determination. 



3.2. Orbital solution 

Using all the available RVs of the secondary component, we 
computed an SB 1 orbital solution. For this purpose, we slightly 
shifted the RVs evaluated on the REOSC spectra as well as the 
data taken from Fullerton ( 1990 ) to account for a slight shift in 
the systemic velocities of the Civ ^ 5801, 5812, Hei A 5876 
lines used for the REOSC and Fullerton ( 1990) data on the one 
hand and the set of lines used to evaluate the secondary RVs from 
our new spectra on the other hand. Indeed, the six lines that we 
use in our determination of the secondary's RV yield systemic 
velocities that differ in the most extreme case by 16 km s _1 (C m 
,15696 vs. Civ ,15801). Since the REOSC and Fullerton ( fT990j ) 
RVs are drawn from a subset of the six lines, we have corrected 
these RVs by +4.3 km s , to avoid biasing our orbital solution. 

The orbital solution was computed using the Liege Orbital 
Solution Package (LOSP) cod^](Sana et al. 120061 ) which is an 
improved version of the code originally proposed by Wolfe et al. 
(fT967l 

For the SB 1 orbital solution, the best-fit orbital period was 
found to be 3163 + 12 days. This solution yields an eccentric- 
ity of 0.69. We next used all 91 available measurements of 
the RVs of both stars (secondary and primary), to compute an 
SB2 orbital solution (see Fig.|2]). The best orbital period is now 
3146 + 18 days and the eccentricity is essentially unchanged 
(e = 0.70). The parameters inferred from this orbital solution 
are listed in TableQ] 



Table 1. Orbital solution of 9 Sgr 



All SB2 data 


Period (days) 


3145.9 ± 17.6 


T (HJD-2400000) 


46927 ± 35 


e 


0.695 ±0.011 


0>C) 


26.4 ± 1.3 


AT, (kms" 1 ) 


25.7 ±0.6 


K 2 (kms" 1 ) 


38.8 ±0.9 


7i (kms -1 ) 


15.8 ±0.4 


y 2 (kms" 1 ) 


5.6 ±0.4 


a sin i (R e ) 


2874 ± 83 


q = OT2/ m i 


0.66 ± 0.02 


mi sin 3 (M e ) 


19.5 ± 1.4 


m 2 sin 3 i (M e ) 


12.9 ±0.9 


rms(0-C)i (kms" 1 ) 


4.2 


rms(0-C) 2 (kms" 1 ) 


3.2 



Notes. To stands for the time of periastron passage, oj is the primary 
star's longitude of periastron measured from the ascending node of the 
orbit. The quoted uncertainties correspond to 1-cr. 



3 This designation might seem strange at first sight, but it is in line 
with our orbital solution (see Sect. |3.2fr . In 9 Sgr, the star with the 
strongest signature in Hei is the less massive and cooler component 
of the system. 



4 The uncertainty on the orbital frequency was estimated as one tenth 
of the natural width of the peaks in the Fourier periodogram as com- 
puted for the full set of secondary RVs. 

5 http://stafF.science.uva.nl/~hsana/losp.html 
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Fig. 2. Radial velocity curve of 9 Sgr. Open and filled symbols 
stand for the RVs of the primary and secondary star respectively. 
The different symbols indicate data taken with different instru- 
ments. The secondary RVs from Fullerton (1990 ) were not in- 
cluded in the computation of the orbital solution. The top panel 
yields the RV curve folded in orbital phase, whilst the lower 
panel illustrates the distribution as a function of time (given in 
years on the top axis and HJD - 2400000 on the bottom axis) of 
all RV measurements used in this paper. 



We note the large difference in systemic velocity between 
the two stars. This is a common feature in massive binaries (e.g. 
Rauw et al. 2002b ) and is usually interpreted as a signature of the 
stellar wind velocity fields affecting the positions of the lines. 
In the present case, this shift can also, at least partially, arise 
from the fact that we use different sets of lines to infer RVi and 
RV2. Indeed, as pointed out above, there are systematic differ- 
ences between the RVs of some of the lines. These differences 
remain roughly constant over the orbital cycle, except for the 
primary's Siiv 41 16 emission and Nv A 4603 absorption lines 
near periastron passage. In fact, at these phases, the Si iv line ap- 
pears red-shifted by +14km s _1 compared to RV 1 whilst theNv 



line is simultaneously blue-shifted by -1 1 km s _1 again with re- 
spect to RVi. This situation could indicate a perturbation of the 
primary's stellar wind due to enhanced wind-wind interactions 
when the stars are at periastron. The immediate consequence of 
this shift is the increase of the cr\ dispersion around the phases 
of periastron passage (see Table A.l I. We do not expect these ef- 



fects to have a major impact on either the amplitudes of the RV 
curves or the results of the disentangling (see Sect. 3.3 1. 

Based on our ephemerides, we predict that the next perias- 
tron passage of 9 Sgr should occur in mid March 2013. 

3.3. Disentangling 

We used our disentangling code based on the method of 
Gonzalez & Levato ( 120061 ) to separate the spectra of the sec- 
ondary and primary component. This method iteratively recon- 
structs the secondary and primary spectra and allows in prin- 
ciple to simultaneously recover the RVs of both components. 
However, owing to its rather modest maximum RV separation, 
9 Sgr is definitely a very difficult case for any disentangling 
method. Indeed, in this system, the spectral lines are clearly 
double only at phases around the secondary's descending node 
(i.e. around periastron passage). The opposite configuration (sec- 
ondary lines red-shifted) does not produce a sufficient RV excur- 
sion to separate the cores of the lines. This situation is expected 
to lead to an imperfect spectral disentangling, especially for the 
H 1 and He 11 lines which have broad wings and do not show clear 
SB2 signatures around the secondary's descending node either 
(see Fig.[T]i. Under these circumstances, it is not obvious that 
the disentangling method can actually improve the RV measure- 
ments. 

To deal with as homogeneous a dataset as possible (in terms 
of spectral resolution and wavelength coverage) that covers most 
orbital phases, we applied the disentangling algorithm to our set 
of FEROS spectra only. We attempted disentangling with differ- 
ent treatments of the RVs, either keeping them fixed for one or 
both stars over part of the iterative process, or allowing the code 
to modify them. The result of our numerous tests is that, in the 
specific case of 9 Sgr, disentangling does not improve the RVs 
and the best results on the spectral reconstruction are actually 
achieved by adopting the RVs determined from direct measure- 
ments of the lines (see above). We then applied the disentangling 
code on four wavelength domains: 4000 - 4360, 4450 - 4750, 
5170 - 5470, and 5580 - 5890A. The results for two regions are 
illustrated in Fig. [3] The disentangled spectra in Fig.[3]were not 
corrected for the optical brightness ratio and their normalization 
hence refers to the combined continua of both stars. 



3.4. Spectral classification and brightness ratio 

Unfortunately, the disentangling does not allow us to well re- 
cover the Hen lines in the spectra of the components of 9 Sgr. 
This is a serious handicap for spectral classification of O-type 
stars which relies to a large extent on the relative strength of 
He 1 and Hen lines. Nonetheless, Fig.[3]reveals a number of in- 
teresting features of the components of 9 Sgr. For instance, we 
observe weak He 1 absorptions in the spectrum of the primary, 
whilst these lines are much stronger in the spectrum of the sec- 
ondary. Despite the uncertainties that affect the Hen lines, we 
can obtain some useful information from these lines as well. 
First, He 11 A 4686 seems to be in relatively strong absorption in 
both stars, suggesting a main-sequence luminosity class for both 
of them. Second, Hen A 4542 appears stronger than Hei A 4471 
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Fig. 3. Disentangled spectra of the secondary and primary com- 
ponents of 9 Sgr in the 4450 - 4750 and 5580 - 5890 A regions. 
The normalization of the spectra is with respect to the level of 
the observed spectra and the primary spectrum is shifted verti- 
cally by 0.1 continuum units. Note that we excluded some parts 
of the yellow spectrum that are affected by strong diffuse inter- 
stellar bands. Note also that, because of the low RV amplitudes, 
there are a number of artefacts, such as the wings of strong and 
broad lines that are not well represented. 



in both stars, indicating an early spectral type for both of them 
(Walborn & Fitzpatrick [T990l . 

The N m AA 4634 - 40 lines are apparently in emission in 
both stars, indicating that both should have an ((f)) tag. The ab- 
sorption lines due to the highest ionization stages of nitrogen 
(N iv and N v) which are absent or very weak in the secondary's 
spectrum, are rather strong in the primary spectrum. Therefore, 
the primary is likely the intrinsically hotter component of the 
system. In this respect, we note also the presence of moder- 
ately strong Si iv /1/14088, 4116 emissions and of a very weak 
Niv A 4058 emission in the spectrum of the primary. Overall, 
the primary spectrum resembles that of HD 93128 classified as 
03.5V((f + )) in Walborn et al. d2002b and we thus adopt the 
same classification for the primary star. The moderately strong 
C in A 5696 emission in the secondary spectrum and the rather 
strong O m and C iv absorption lines are indicative of a spectral 
type later than 04 (Walborn 1980). Comparing the overall as- 
pect of the secondary spectrum with the catalogue of Walborn 
et al. ( 120021 1. we find that the best agreement is found with the 
spectrum of HD 46150, an 05 V((f)) star according to Walborn 



et al. (120021) reclassified to 05.5 V((f)) by Mahy et al. (2009), 
and we thus adopt an 05-5.5V((f)) spectral classification for 
the secondary. Note that, according to the spectroscopic masses 
quoted by Martins et al. ( 120051 ). the dynamical mass ratio of 
q = 0.66 + 0.02 inferred from our orbital solution, is in rea- 
sonable agreement with the expected mass ratio (0.66 - 0.73) 
between an 05-5.5V secondary and an 03.5 V primary. Also, 
comparing our minimum masses from TablefTjwith the spectro- 
scopic masses quoted by Martins et al. (120051 ). we estimate an 
orbital inclination of 45° ± 1°. Given the rather wide separation 
between the two components of the binary, interferometry can 
resolve the system and should provide a determination of the 
separation of the stars projected on the sky. Indeed, assuming 
i = 45° and adopting a distance of 1.79kpc (see below), we 
predict angular separations of 3 mas at periastron and 17 mas at 
apastron. These numbers are within reach of modern interfero- 
metric facilities such as the VLTI. Combined with our orbital 
solution and an estimate of the distance, such a measurement 
yields, in principle, an independent determination of the orbital 
inclination and hence the absolute masses of the binary compo- 
nents. 

From a preliminary attempt to disentangle the secondary and 
primary spectra, Rauw et al. ( 12005 1 inferred spectral types 04 
for the primary and 07-8 for the secondary. Our new results are 
quite different. The reason for this discrepancy stems from the 
fact that our previous analysis was restricted to the region of 
the He i A 5876 and C iv /LI 5801, 5812 lines and we erroneously 
assumed the star with the strongest signature in these lines to be 
the hotter component of the system. As we have seen above, this 
is not the case. 

We measured the equivalent widths (EWs) of some He i and 
some metallic lines on the disentangled spectra of the two com- 
ponents (but referring to the combined continuum of the two 
stars). Comparing these numbers with the EWs observed in typ- 
ical 05 and 05.5 stars as tabulated by Conti (11974b and Conti 
& Frost (11 9771 ), we estimated the dilution of the secondary's 
spectral signature by the light of the primary. Assuming an 05 
spectral type for the secondary, we find that this star contributes 
0.50 + 0.05 of the total light in the optical domain. If we assume 
an 05.5 type instead, the contribution would be 0.39 ± 0.04. In 
other words, the secondary star should roughly contribute be- 
tween 40 and 50% of the total light. This would then imply an 
optical brightness ratio of about 1.0 - 1.5 for the primary com- 
pared to the secondar^ This result is in reasonable agreement 
with the expected difference in M v between an 03.5 V and an 
05-5.5 V star. Indeed, from Table 4 in Martins et al. ( 120051 ). we 
expect an optical brightness ratio in the range 1.5 - 1.7. 

Unfortunately, there are not enough EWs of the earliest O- 
type stars tabulated in the literature, so that we cannot perform 
an independent evaluation of the brightness ratio on the EWs of 
the primary star. 

If we assume that both stars have typical absolute magni- 
tudes for their spectral types, we estimate a combined My of 
-6.20. With an apparent V magnitude of 5.93 and B - V — 0.0, 
we then estimate a distance modulus of 1 1 .26 (d = 1 .79 kpc), in 
extremely good agreement with the distance of the open cluster 
NGC 6530 (1.78 + 0.08 kpc) as inferred by Sung et al. d2000l) . 



6 In our preliminary analysis (Rauw et al. 2005 ), we inferred a bright- 
ness ratio of 4 ± 1 based on the equivalent width ratio of the He i A 5876 
and Civ /M5801, 5812 lines. The discrepancy with our current result 
is entirely due to the different spectral types assumed in our previous 
work (see above). 
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We used the disentangled spectra to estimate the projected 
rotational velocities v eq sin i of the stars. For this purpose, we ap- 
plied the Fourier transform method (Gray 2005 Simon-Diaz & 
Herrero|2007]) to the O m A 5592, C rv AA5801, 5812 lines of the 
disentangled spectra. Indeed, these lines are the ones that are best 
reconstructed by the disentangling method. In some cases, the 
Fourier transforms of the lines are not very well represented by 
the theoretical Fourier transform of a rotational broadening func- 
tion. However, by concentrating on the position of the first 'zero' 
of the Fourier transform, we estimate v eq sin i = (87 ±16) km s 
for the primary and (57+5) km s _1 for the secondary. The quoted 
errors correspond to the dispersions about the mean of the values 
determined for the three lines and does not account for system- 
atic uncertainties due e.g. to the presence of macroturbulence 
that could impact the measurement of the rotational velocity. 

As a consistency check of the above results, we have simu- 
lated the synthetic spectra of an 03.5 V star and an 05 V star. 
For this purpose we have used the CMFGEN non-LTE atmo- 
sphere model (Hillier & Miller 1998) which accounts for the 
presence of a stellar wind and includes the effects of line blan- 
keting. In this model atmosphere code, the equations of radiative 
transfer are treated in the comoving frame and the hydrodynamic 
structure of the atmosphere is built from the TLUSTY models 
of the OSTAR2002 grid of Lanz & Hubeny d2003l >. We caution 
that our models are by no means adjusted to our disentangled 
spectra. Indeed, as pointed out above, there are artefacts in some 
lines of the disentangled spectra and furthermore, in our spe- 
cific case, we lack reliable information on the mass-loss rates of 
the stars. However, these models were computed adopting some 
mean parameters of 03.5 and 05 dwarfs taken from Martins et 
al. ( 120051 ). For the primary, we assume T e $ = 43850 K, logg = 
3.92 (in cgs units), M = 9 x IV 1 M^yr 1 , v M = 3500 kms" 1 , 
whilst the secondary parameters were taken as r e ff = 40850 K, 



logg = 3.92, M = 5 x l(r'M G yr- 



3 100 km s" 1 . Each 



spectrum was further broadened by the observationally deter- 
mined v sin/. These synthetic spectra allowed us to check the 
brightness ratio and the behaviour of the He n lines. Concerning 
the first point, we evaluated the dilution factors of the spectra 
of each component by comparing its observed EWs to the simu- 
lated ones. In this way, we infer an optical brightness ratio (pri- 
mary/secondary) of 0.87 ±0.40 in reasonable agreement with the 
values derived hereabove. As to the second issue, by combining 
the two synthetic spectra shifted by the RVs at maximum sepa- 
ration, we recover the observed behaviours described in Sect. [3] 
and Fig.[T[ the Hei, Om and Civ lines are indeed found to dis- 
play an SB2 signature whilst the He n lines, which have roughly 
equal strength in both stars, do not display a clear binary signa- 
ture. Finally, we note that, in agreement with our observations, 
Cm A 5696 emission is significantly present only in the synthetic 
05 V spectrum, whilst the CMFGEN spectra do not reproduce 
the observed N v absorptions in the primary spectrum. Overall, 
the agreement is quite reasonable and lends further support to 
our spectral classification and optical brightness ratio near unity. 

4. Discussion 

Wind-wind interactions in massive binaries can leave their sig- 
nature over a wide range of wavelength (see e.g. Rauw 120101 
for a review). Among the most common signatures of this phe- 
nomenon are optical emission line variations, (variable) excess 
X-ray emission and non-thermal radio emission. 

In close O-type binaries, the He n A 4686 and Ha emission 
lines are frequently seen to display phase-locked variations due 
to the wind-wind collision (e.g. Linder et al. 2008 and references 



therein). However, in 9 Sgr, these lines remain in absorption all 
over the orbital cycle and display little profile variations that 
are most easily explained by a partial deblending of the lines 
near periastron. In the Ha line, the only emission is a weak nar- 
row nebular emission arising from the Lagoon Nebula in which 
9 Sgr is embedded. This lack of optical emission lines is most 
likely due to the fact that the wind interaction in 9 Sgr remains 
in the adiabatic regime over the entire orbital cycle. Indeed, in 
this case, the post-shock gas would not undergo radiative cool- 
ing and would not produce a high density region of sufficiently 
low temperature to produce significant line emission. There are 
two reasons why the wind interaction in 9 Sgr would not be in 
the radiative regime. First the components of 9 Sgr are still on the 
main-sequence and therefore feature rather low mass-loss rates. 
Second the separation between the two stars is quite large, even 
at periastron, where r sin/ = 877 7? Q = 4AU, r being the in- 
stantaneous separation between the stars. This implies that the 
pre-shock wind density is likely very low, preventing radiative 
cooling from being effective. 

The XMM-Newton observation discussed by Rauw et al. 
(2002a} was obtained on March 8, 2001, i.e. at orbital phase 
0.605. We re-analysed this dataset in the light of the results 
presented above. The EPIC spectra can be fitted with a two- 
temperature thermal plasma model featuring a soft component 
with kT\ = 0.28 keV and a hard component at kT2 = 2.5 keV. 
The best-fit model has an X-ray flux (corrected for the interstel- 
lar absorption) of 4.3 x 10~ 12 erg cm" 2 s over the 0.5 - lOkeV 
spectral domain. This number translates into an X-ray luminos- 
ity of 1.65 x 10 33 ergs _1 assuming a distance of 1.79 kpc. As 
pointed out by Rauw et al. (2002a), this X-ray luminosity is not 
exceptionally large, and does not provide direct evidence for a 
contribution of the wind-wind interaction. However, this obser- 
vation was taken at phase 0.605, i.e. near apastron. If we as- 
sume that the wind interaction region in 9 Sgr is in the adiabatic 
regime (see above), then the X-ray emission due to the colliding 
winds should vary as l/r (Stevens et al. 11992b . Such a variation 
is indeed observed in some wide O-type binaries, a spectacu- 
lar example being the 2.35 yr period system Cyg OB2 #9 (Naze 
et al. 2012, submitted). If the entire X-ray emission of 9 Sgr is 
due to the colliding wind interaction, a l/r scaling would im- 
ply that the X-ray flux of 9 Sgr should increase by a factor 5.4 
near periastron (r = 0.305 a) compared to our XMM-Newton ob- 
servation (r = 1 .64 a). Alternatively, if we assume that only the 
harder spectral component seen in our EPIC spectra arises from 
the wind- wind interaction and hence undergoes the l/r varia- 
tion, we would still expect an increase of the observed X-ray 
flux by 50% near periastron. A dedicated monitoring of the X- 
ray emission of the system around its next periastron would be 
most helpful to clarify this issue. 

As to the presence of relativistic electrons accelerated in the 
wind-wind collision zone, we stress that a VLA radio observa- 
tion taken simultaneously with our XMM observation revealed 
a negative spectral index typical of non-thermal radio emission 
(Rauw et al.|2002a). Therefore, the shock between the winds is 
capable of accelerating electrons to relativistic velocities, even 
at phases close to apastron when the separation is largest. On the 
basis of the flux density measurements and on the radio spectral 
index between 3.6 and 6 cm determined by Rauw et al. (f2002a), 
we estimate that the synchrotron luminosity at the epoch of the 
corresponding radio observations is of the order of 10 29 ergs _1 . 
On the other hand, the cumulated kinetic power of the two stellar 
winds calculated from the quantities given at the end of Sect. 3.4 



should be of the order of 5 x 10 ergs . This means that, at 
the orbital phase of the radio observation, a fraction of the or- 
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der of 2 X 10 of the kinetic power was converted into syn- 
chrotron radio luminosity. Considering first that at most a few 
percent of the kinetic power are injected in the wind-wind inter- 
action, and second that a similar fraction of the injected power 
is used for particle acceleration, a power of about 10 33 erg s _1 is 
available for relativistic particles, among which some are likely 
involved in non-thermal radiation processes. The measured syn- 
chrotron luminosity corresponds to the fractional power initially 
injected into relativistic electrons that is not radiated through 
inverse Compton scattering (the dominant energy loss process 
for relativistic electrons in colliding- wind binaries; see Pittard & 
Dougherty 2006, De Becker 2007 ). Assuming that the electron- 
to-proton ratio among the population of relativistic particles is 
of the order of 0.01 (as in the case of Galactic cosmic rays), con- 
verting about one percent of the energy of the relativistic elec- 
trons into synchrotron radio emission would be sufficient to ac- 
count for the observations. Such a ratio is in agreement with the 
expected magnetic to radiation field energy densities in typical 
colliding-wind binaries (Pittard & Dougherty 12006b . Therefore, 
within our current understanding of non-thermal processes in 
massive binaries, the above considerations show that the syn- 
chrotron radiation of 9 Sgr could indeed be entirely provided via 
an acceleration process taking place in the wind-wind interaction 
region of this binary system. The current uncertainties on the de- 
tails of some of the energy conversion processes involved in the 
acceleration processes prevent us from doing a more extensive 
discussion in the present study. 

5. Conclusions 

We have clearly established the multiplicity of 9 Sgr and de- 
rived its first ever orbital solution. With an orbital period of 
8.6 years, 9 Sgr is one of a few O + O spectroscopic bina- 
ries with orbital periods longer than one year. Other well estab- 
lished O-type spectroscopic binaries in this part of the parame- 
ter space are HD 15558 (P = 440 days, De Becker et al.[2006), 
Cyg OB2 #9 (P =858 days, Naze et al. l2010al 2012) and 15 Mon 
(P = 9247 days, Gies et al. 1 1 9931 >. It is very likely that a number 
of similar systems remain to be discovered in this observational 
no-man's land. As demonstrated by our present study, exploring 
this part of the parameter space requires high-resolution spectro- 
scopic monitoring over long periods of time. 

Our results support the current paradigm that attributes the 
presence of synchrotron radio emission in early-type stars to a 
wind-wind collision in a binary system. It is striking though that 
O-type non-thermal radio emitters discovered to be binaries over 
the last decade span a wide range of parameters. Indeed, the list 
of such stars includes systems with surprisingly short orbital pe- 
riods such as CygOB2 #8a (P = 21.9 days, De Becker et al. 
2004) as well as much wider systems, the most extreme one so 
far being 9 Sgr. At the same time, the list of non-thermal radio 
O-type stars includes main-sequence stars as well as giants or 
supergiants, i.e. stars having very different mass-loss rates. A 
characterization of these systems, both in terms of orbital and 
stellar parameters, is mandatory to properly understand the pro- 
duction of non-thermal electrons in O-star binaries and quantify 
their contribution to the acceleration of cosmic particles. 
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Table A.l. Journal of our observations of 9 Sgr 



HJD 


Exp. time 

(s) 


Instrument 


(kms-») 


(kms- 1 ) 


RV 2 
(kms- 1 ) 


°"2 

(km s~ 


2451327.93054 


450 


F 


7.1 


5.4 


18.7 


2.9 


2451669.71311 


420 


F 


11.5 


3.7 


14.5 


1.8 


2451670.69770 


240 


F 


12.3 


2.6 


16.0 


3.6 


2451671.71623 


360 


F 


13.4 


4.6 


14.3 


3.4 


2451672.68625 


300 


F 


12.5 


3.1 


14.3 


2.7 


2451673.68442 


480 


F 


6.0 


16.2 


11.8 


3.4 


2452037.72550 


600 


F 


8.4 


8.5 


11.0 


4.9 


2452039.92845 


300 


F 


21.6 


2.0 


10.0 


4.0 


2452040.68815 


540 


F 


18.6 


6.5 


12.0 


2.8 


2452102.71686 


90 


U 


8.6 


0.7 


5.1 


- 


2452335.89786 


300 


F 


15.3 


9.7 


7.9 


2.2 


2452336.88868 


600 


F 


12.9 


4.8 


7.8 


8.9 


2452337.89775 


600 


F 


17.3 


1.7 


8.5 


4.7 


2452338.90100 


600 


F 


13.1 


2.8 


9.3 


3.2 


2452339.91296 


600 


F 


11.5 


4.5 


7.0 


3.8 


2452353.88979 


90 


Em 


19.5 


17.2 


1.5 


4.1 


2452381.87295 


1200 


F 


14.8 


0.8 


6.3 


3.7 


2452745.86397 


450 


Es 


19.4 


2.7 


-3.7 


9.5 


2452748.88352 


540 


Es 


18.0 


5.0 


-8.9 


7.8 


2452782.76326 


700 


F 


23.1 


3.1 


-8.6 


5.7 


2452783.92065 


900 


F 


23.1 


4.0 


-4.0 


6.4 


2452784.93547 


600 


F 


29.3 


6.6 


-5.7 


5.6 


2453130.93332 


240 


F 


54.9 


13.3 


-52.0 


7.1 


2453133.88227 


240 


F 


57.9 


11.1 


-51.0 


6.0 


2453139.68916 


450 


F 


51.5 


15.1 


-55.0 


4.7 


2453159.73683 


450 


F 


54.6 


17.7 


-52.1 


5.2 


2453196.58561 


450 


F 


56.5 


17.4 


-57.2 


4.6 


2453226.64288 


900 


F 


55.9 


23.8 


-51.0 


7.4 


2453231.59984 


450 


F 


48.4 


17.1 


-49.1 


8.6 


2453240.59441 


450 


F 


47.5 


13.4 


-44.1 


6.6 


2453243.52397 


450 


F 


48.9 


9.9 


-42.1 


6.4 


2453246.58287 


450 


F 


47.2 


14.3 


-42.3 


7.3 


2453255.52210 


450 


F 


43.0 


10.6 


-38.4 


7.6 


2453271.52285 


300 


F 


42.4 


11.3 


-34.0 


6.9 


2453273.51418 


450 


F 


40.0 


9.6 


-37.0 


7.2 


2453288.60031 


840 


Es 


34.2 


9.4 


-27.6 


11.2 


2453289.59542 


600 


Es 


28.9 


8.3 


-34.1 


36.6 


2453503.84288 


450 


F 


8.6 


2.1 


14.6 


3.1 


2453509.88418 


300 


F 


10.5 


1.9 


14.7 


2.8 


2453521.74907 


450 


F 


10.9 


4.3 


12.9 


3.3 


2453561.74156 


450 


F 


11.2 


1.9 


15.8 


3.0 


2453597.61875 


450 


F 


8.1 


0.7 


18.2 


2.7 


2453664.52625 


450 


F 


7.1 


5.8 


19.6 


3.2 


2453681.54040 


450 


F 


6.3 


2.1 


18.8 


3.7 


2453796.89444 


300 


F 


5.5 


4.8 


18.0 


2.3 


2453855.94364 


300 


F 


4.0 


4.0 


22.5 


4.6 


2453861.92973 


300 


F 


6.9 


4.7 


20.6 


4.7 


2453874.85744 


600 


F 


5.2 


4.0 


21.6 


3.4 


2453910.64138 


800 


F 


5.7 


1.3 


20.3 


2.4 


2453919.56510 


600 


F 


7.5 


1.4 


23.7 


8.6 


2453958.47874 


600 


F 


3.6 


5.6 


20.7 


4.2 


2454013.52591 


414 


F 


4.1 


2.8 


24.2 


4.5 


2454184.88276 


552 


F 


6.7 


4.0 


17.7 


2.8 


2454193.97917 


900 


Es 


-3.4 


25.2 


15.4 


5.2 


2454210.93497 


120 


F 


4.5 


3.0 


22.1 


6.2 


2454212.82582 


138 


F 


4.8 


2.6 


19.6 


3.6 


2454267.72914 


138 


F 


9.4 


2.9 


18.2 


4.9 


2454274.86447 


138 


F 


5.4 


1.8 


17.3 


3.1 



Notes. The codes for the instrument correspond to Coralie (C), FEROS (F), EMMI (Em), Espresso (Es) and UVES (U). The columns labelled cr { 
and <T2 provide the 1-cr dispersion about the mean RV of the lines that are used to evaluate the radial velocity of a given component (see Sect.[3]l. 
For a few UVES observations, no number is provided because the instrumental set-up covered only one of the lines used here. 



8 



Rauw et al.: The O-type binary 9 Sgr 

Table A.l. continued 



HJD 


Exp. time 


Instrument 


RVi 




RV 2 


0~2 




(s) 




(kms ') 


(kms ') 


(kms ') 


(kms ') 


2454333.53105 


138 


F 


4.3 


2.6 


19.1 


2.2 


2454372.53271 


138 


F 


8.2 


5.9 


19.9 


1.9 


2454378.49630 


150 


U 


13.7 


6.1 


18.4 


8.3 


2454540.82656 


150 


U 


12.2 


4.8 


17.6 


4.7 


2454595.83785 


318 


c 


13.6 


15.4 


17.3 


3.4 


2454658.67968 


318 


c 


-5.8 


7.8 


14.4 


2.7 


2454725.48261 


565 


c 


7.1 


3.3 


15.9 


2.1 


2454765.52334 


150 


u 


10.2 


2.7 


19.5 


6.0 


2454781.50546 


565 


c 


10.7 


3.0 


15.0 


4.5 


2454874.88694 


600 


c 


3.1 


15.1 


14.5 


3.0 


2454905.81592 


150 


u 


16.5 


8.2 


17.5 


2.8 


2454935.78120 


624 


c 


6.4 


9.4 


14.8 


3.2 


2454941.90945 


30 


u 


12.4 


5.1 


7.8 


- 


2454956.93860 


250 


F 


7.4 


6.6 


13.0 


2.9 


2454995.62614 


600 


c 


10.0 


11.5 


13.0 


2.4 


2455064.56491 


624 


c 


15.8 


4.6 


12.2 


3.8 


2455132.49691 


600 


c 


16.4 


4.8 


12.2 


2.5 


2455247.89285 


660 


c 


6.8 


7.8 


11.6 


1.8 


2455284.89819 


30 


u 


11.5 


3.3 


7.3 


— 


2455311.89860 


900 


c 


14.8 


0.9 


11.8 


2.9 


2455365.75019 


624 


L 


7.1 


n a 

1 A 


9.8 


1 A 

1.4 


2455443.58301 


624 


C 


19.3 


10.4 


8.8 


3.3 


2455478.53930 


30 


U 


10.9 


1.7 


17.2 




2455487.49425 


600 


C 


15.1 


6.9 


7.9 


4.8 


2455622.83675 


30 


U 


8.7 


2.2 


16.5 




2455656.88117 


624 


c 


12.1 


3.6 


7.8 


2.8 


2455668.90704 


60 


u 


11.0 


3.9 


-1.4 




2455724.71507 


1200 


Es 


30.7 


10.2 


13.2 


7.7 



